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. Abstract 

! 

An upper bound on the ergodic capacity of MIMO channels was introduced recently 
in [I] . This upper bound amounts to the maximization on the simplex of some multilinear 
polynomial p(Ai, A n ) with non-negative coefficients. In general, such maximizations prob- 
lems are NP-HARD. But if say, the functional log(p) is concave on the simplex and can be 
| efficiently evaluated, then the maximization can also be done efficiently. Such log-concavity 

was conjectured in [TJ. We give in this paper self-contained proof of the conjecture, based 
on the theory of H-Stable polynomials. 
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1 The conjecture 

O . 

Let B be M x M matrix. Recall the definition of the permanent : 

Per{B)= £ II ^(*,<7(i)). 
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<TgS M l<i<M 

The following Conjecture was posed in pp. 

Conjecture 1.1: Let AkMxJV,M<JV matrix with non-negative entries. We denote as 
As a submatrix 

A s = {A(i,j) :l<i<m;jeSc{l, N}. 
Define the following multi-linear polynomial with non-negative coefficients 

F A (X 1 ,...,X N )= ]T Per(A s )Y[X 3 . (1) 

\S\=M,SC{1,...,N} jes 

Then the functional log(FA) is concave on = {(Ai, Ajv) : Aj > 0, 1 < j < N}. | 

We present in this paper a proof of Conjecture (1 1.11) . Actually we prove that the polynomial 
Fa is either zero or H-Stable. 
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2 H-Stable polynomials 



To make this note self-contained, we present in this section proofs of a few necessary results. 
The reader may consult J5^ and J3}/ for the further reading and references. 

We denote as Hom+(m,n) a convex closed cone of homogeneous polynomials with non- 
negative coefficients of degree n in m variables and with non-negative coefficients; as R™ a 
convex closed cone of non-negative vectors in R m and as a convex open cone of positive 
vectors in R m . 

Definition 2.1: A homogeneous polynomial p G Hom+(m,n) is called H-Stable if 

\p(z 1 , ...,z m )\ > 0;Re(zi) > 0, 1 < i < m; 

is called H-SStable if \p{z\, ...,z m )\ > provided that 

Re(zi) > 0, 1 < i < m and < J2i<m<m R e ( z i)- 

I 

Example 2.2: Consider a bivariate homogeneous polynomial p G Hom + (2,n), p{z\,Z2) = 
(z2) n -P(f^)i where P is some univariate polynomial. Then p is H-Stable iff the roots of P are 
non-positive real numbers. This assertion is just a rephrasing of the next set equality: 

C - {— : Re(z 1 ),Re(z 2 ) > 0} = {x G R : x < 0}. 

In other words 

P(t)=a Yl (t + ai) : ai> 0,1 <i < k;a > 0. 

l<i<k<n 

Which gives the following expression for the bivariate homogeneous polynomial p: 

l<i<k<n 

I 

Fact 2.3: Let p G Hom+(m,n) be H-Stable. Then log(p) is concave on R+. 

Proof: Consider two vectors X, Y G .R+ such that their sum X + Y G R™ has all positive 
coordinates. It is sufficient to prove that the bivariate homogeneous polynomial q G Hom+(2, n) 

q{t,s)=p{tX + sY), 

is log-concave on Clearly, the polynomial q is H-Stable. Therefore, using Example (|2.2j) . 
we get that 

log(g(t, s)) = log(a) + (n — k) log(s) + V] log(t + ajs) : > 0, 1 < i < A;; a > 0. 

l<i<fc<n 

The log-concavity of q follows now from the concavity of the logarithm on [0, oo). | 
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Remark 2.4: Since the polynomial p is homogeneous of degree n hence, by the standard 
argument, the function p» is concave on R™ as well. | 

Fact 2.5: Let p £ Hom + (m,n) be H-Stable and x^ > 0, 1 < i < m then the following 
inequality holds 

\p(xi + iyi, ...,x m + iy m )\ > p(x±, ...,x m ) (2) 

Proof: Consider without loss of generality the positive case Xi > 0, 1 < i < m. Then there 
exists a positive real number \i > such that + \iX{ > 0, 1 < i < m. It follows from 
Example ()2. 2 p that for all complex numbers z € C 

p{zx x + (yi +fj,xx), ...,x m + z(y m + fj,x m ) = p(x u ...,x m ) J[ {z + cii);ai > 0, 1 < i < m. 

\<i<n 

Thus 

p(zxi+yi,...,zx m + y m ) =p(xi,...,x m ) J| (z + ai - p) 

l<i<n 

We get, using the homogeniuty of the polynomial p, that 

p(x 1 + iyi,...,x m + iy m ) = p(xi,...,x m ) J| (1 + i(aj - //)). 

i<i<« 

As | ]li<j<n(l + K a j ~ A 4 )) I ^ 1 t n i s proves that the inequality ([2]) holds. I 

Corollary 2.6: A nonzero polynomial p £ Hom + (m,n) is H-Stable if and only the inequlity 
(3) holds. 

Corollary 2.7: Let pi £ Hom + (m,n) be a sequence of H-Stable polynomials and p = 
linij^oop,;. Then p is either zero or H-Stable. 

Some readers might recognize Corollary \2. 7| j as a particular case of A. Hurwitz's theorem on 
limits of sequences of nowhere zero analytical functions. Our proof below is elementary. 

Proof: Suppose that p is not zero. Since p £ Hom + (m, n) hence p(x\, . . . , x m ) > if Xj > : 
1 < j < m. As the polynomials pi are H-Stable therefore |pj(Z)| > \pi{Re(Z)) \ : Re(Z) £ 
R^ + . Taking the limits we get that \p(Z)\ > \p(Re(Z)) | > : Re(Z) £ R% + , which means 
that p is H-Stable. | 

We need the following simple yet crucial result. 

Proposition 2.8: Let p € Hom + (m,n) be H-Stable. Then the polynomial pi±\ € Hom + (m — 
l,n-l), 

d 

P(i)(x 2 ,-,x m ) =: —p(0,x 2 ,...,x m ), 

is either zero or H-Stable. 
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Proof: Fix complex numbers z£, 2 < i < m and define the following univariate polynomial 

R(t) = p(t,z 2 ,...,z m ). 
It follows that R'(0) = P(i)( z 2, •••) z m ). We consider two cases. 

First case: the polynomial p £ Hom + (m,n) is H-SStable. In this case the polynomial G 
Hom + (m — 1, n — 1) is H-SStable as well. Indeed, in this case if the real parts Re(zi) > 0, 2 < 
i < m and X^2<i<m R^( z i) > then all the roots ui, ...,u ri _i of the univariate polynomial -R 
have strictly negative real parts: 

R(t) = h Yl {t-Vi),0^heC. 

2<«<n-l 

Therefore 

P(1) (z 2 ,...,z m ) = J R'(o) = / i (-ir- 2 ( n E Kr 1 )^ 

2<i<n-l 2<i<n-l 

as the real part 

**( E N- 1 )= E W>°- 

2<i<n-l 2<i<n-l ' l| 

Second case: the polynomial p £ Hom + (m,n) is H-Stable but not H-SStable. We need to 
approximate p by a sequence of H-SStable polynomials. Here is one natural approach: let A 
be any m x m matrix with positive entries. Define the following polynomials: 

p I+eA (Z) =:p((I + eA)Z),Z eC m . 

Clearly, the for all e > the polynomials pi+ t A G Hom+(m,n) and are H-SStable. 
It follows that polynomials -^pj +e A(0, X2, x m ) are H-SStable as well. Note that 

d 

llU k -^—Pl+eA{0, Z 2 , Zm) = Pm[Z2, Z m ). 

Using Corollary (|2.7p we get that the polynomial pn\ is either H-Stable or zero. | 



3 Proof of the conjecture 

Proof: We will need a few auxiliary polynomials: 

P(x 1 ,...,x M ;X 1 ,...,X N ) = (Xj + E A (h3) x i)- ( 3 ) 

±<3<N l<i<m 

Clearly, the polynomial P £ Hom+(M + N, N) is H-Stable if the entries of the matrix A are 
non-negative. Applying Proposition (|2.8p inductively, we get that the following polynomial 

d m 

i?(Ai, \ N ) = Qx P(X = 0; Ai, X N ) (4) 
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is either zero or H-Stable as well. It is easy to see that 

R(X 1 ,...,X N ) = J2 Per(A s ) JjAj, (5) 
\S\=M,Sc{i,...,N} jeS 

where S = {1, N} — S is the compliment of the set S. 

Now everything is ready for the punch line: the multilinear homogeneous polynomial, 
defined in (pQ), 

F A (X 1 ,...,X N ) = ( J] X i )R((X 1 )-\...,(X N )- 1 ). (6) 

l<i<N 

Recall that the real part Re{z^ 1 ) = f° r an non-zero complex numbers z £ C. Therefore, 
if the real parts Re(Xi) > 0, 1 < i < n then the same is true for the inverses: 

Re{{XiY x ) > 0,1 < i < n. 

This proves that the polynomial Fa is either zero or H-Stable. The log-concavity follows from 
Fact (H3|) . | 

4 Conclusion 

The reader should not be deceived by the simplicity of our proof: very similar arguments are 
behind the breakthrough results in [5], [4], [6]. The reader is advised to read very nice exposition 
in ®- 

Conjecture (11. ip is actually a very profound question. Had it been asked and properly answered 
in 1960- 70s, then the theory of permanents (and of related things like mixed discriminants and 
mixed volumes [6]) could have been very different now. 

Though the "permanental" part in [T] is fairly standard(the authors essentially rediscovered so 

called Godsil-Gutman Formula [8]) it is quite amazing how naturally the permanent enters the 

story. Switching the expectation and the logarithm can be eventful indeed. 

The log-concavity comes up really handily in the optimizational context of [lj. The thing is 

that maximization on the simplex of J2i<i<j<Nb(i,j)xiXj is NP-COMPLETE even when 

6(i,i)e{0,l},l<l<j<JV. 

Our proof is yet another example on when the best answer to a question posed in the real 
numbers domain lies in the complex numbers domain. Yet, we don't exlude a possibility of a 
direct "monstrous" proof. 

References 

[1] X. Gao, B. Jiang, X. Li, A. B. Gershman and M. R. McKay, M Statistical Eigenmode 
Transmission Over Jointly Correlated MIMO Channels, IEEE TRANSACTIONS ON IN- 
FORMATION THEORY, VOL. 55, NO. 8, AUGUST 2009, 3735-3750. 

[2] A. Schrijver, Counting 1-factors in regular bipartite graphs, Journal of Combinatorial 
Theory, Series B 72 (1998) 122-135. 



5 



[3] M. Laurent and A. Schrijver, On Leonid Gurvits' proof for permanents, 2009, 
http://homepages.cwi.nl/~lex/files/perma5.pdf, to appear in American Mathematical 
Monthly. 

[4] S. Friedland and L. Gurvits, Lower Bounds for Partial Matchings in Regular Bipartite 
Graphs and Applications to the Monomer-Dimer Entropy, Combinatorics, Probability and 
Computing, 2008. 

[5] L. Gurvits, Van der Waer den/ Schrijver- Valiant like conjectures and stable (aka hyperbolic) 
homogeneous polynomials: one theorem for all, Electronic Journal of Combinatorics 15 
(2008). 

[6] L. Gurvits, A polynomial-time algorithm to approximate the mixed volume within a simply 
exponential factor. Discrete Comput. Geom. 41 (2009), no. 4, 533-555. 

[7] L. Gurvits, On multivariate Newton-like inequalities, http://arxiv.org/abs/0812.3687, 

[8] C. Godsil and I. Gutman. On the matching polynomial of a graph. Algebraic Methods in 
Graph Theory, pages 2412013249, 1981. 



6 



